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Project  Objectives: 

The  objective  of  this  research  is  to  study  lipid  bilayer  models  by  molecular  dynamics 
simulations.  The  computer  code,  which  has  been  developed  on  the  Cray  XMP-24  at  the 
Naval  Research  Laboratory  (NRL),  is  based  on  the  Monotonic  Lagrangian  Grid  (MLG) 
near  neighbors  algorithm  developed  by  Boris  and  Lambrakos  of  the  Laboratory  for 
Computational  Physics  and  Fluid  Dynamics  at  NRL  [1,2].  The  MLG  algorithm  had  been 
developed  for  simulating  diatomic  systems  such  as  N2  [3].  A  characteristic  deeply 
embedded  in  the  current  realization  of  the  MLG  algorithm  is  two-dimensional  periodicity. 

The  direction  normal  to  the  periodic  plane  can  be  bounded  by  a  reflecting  wall  or  it  may 
be  free.  A  constraint  algorithm  had  been  developed  and  tested  to  maintain  a  fixed 
distance  between  the  constitutive  atoms  of  the  diatomic  molecules. 

There  were  two  paramount  problems  to  solve  in  order  to  adapt  the  molecular 
dynamics  codes  to  simulations  of  systems  of  polyatomic  molecules  such  as  lipid  bilayers. 
First,  it  was  necessary  to  generalize  the  constraint  algorithm  to  the  polyatomic  case.  This 
work  has  been  described  in  publication  #1,  below.  Second,  it  was  desirable  to  develop  an 
efficient  method  for  calculating  the  intramolecular  angle  dependent  functions  for  potentials 
and  forces  of  torsion  and  angle  bending.  The  angle-dependent  force  algorithm  determines  __ 
the  directions  along  which  forces  are  to  be  applied  to  particles  by  projection  methods, 
rather  than  by  calculating  cross  products.  The  projection  code  used  fewer  expensive 
operations:  the  angle-bending  algorithm  is  approximately  7  times  faster  and  the  torsion  ■- 

algorithm  approximately  40%  faster  than  the -cross-product  based  code.  _  This  work  has  _ 

been  completed,  and  a  manuscript  is  in  preparation  for  submission  to  Siam  J.  Sci.  StaL  _ 

Comp.  _  J(j' 
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A  forty  picosecond  simulation  of  a  butane  system  has  been  conducted  to  test  the 
performance  of  the  constraint  and  torsion  algorithms.  The  Ryckaert-Bellemans  torsional 
and  intermolecular  potential  (Lennard-Jones  with  a  cross  section  of  3.92  A  and  well  depth 
of  72  K)  were  used  [4].  The  energy  of  the  system  versus  time  is  shown  in  figure  1.  The 
initial  angular  distribution  was  chosen  with  a  gaussian  random  distribution  about  the  trans 
(0°)  minimum  of  the  torsional  potentials,  and  the  molecules  were  assigned  random  center 
of  mass  velocities.  The  total  energy  is  reasonably  well  conserved;  improvements  are 
possible  by  using  a  shorter  time  step  or  by  using  multiple  time  steps  to  separate  rapid  and 
slow  degrees  of  freedom.  The  development  of  the  distribution  of  torsion  angles  towards 
the  equilibrium  Boltzmann  distribution  is  shown  in  figure  2  at  0,  1,  5,  10,  25,  and  40  ps. 

Simulations  of  a  C10  bilayer  are  underway. 
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